Quantum computing of delocalization in small-world networks 
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We study a quantum small-world network with disorder and show that the system exhibits a 
delocalization transition. A quantum algorithm is built up which simulates the evolution operator 
of the model in a polynomial number of gates for exponential number of vertices in the network. 
The total computational gain is shown to depend on the parameters of the network and a larger 
than quadratic speed-up can be reached. We also investigate the robustness of the algorithm in 
presence of imperfections. 
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Recently, much attention has been attracted to the 
study of small-world networks Q] . They have been shown 
to describe social and biological networks, Internet con- 
nections, airline flights and other complex networks. In 
such systems, it is possible to go from a given point to 
any other through only a small number of links. Well- 
established classical models have been proposed and ana- 
lyzed by statistical methods. The study of quantum net- 
works with the same property has started only recently, 
showing that these systems present interesting features 
related to quantum transport, delocalization [3, |3j and 
fast diffusion £|. 

In parallel, the development of quantum information 
and computation has become more and more important 
0- In particular, the study of quantum computers has 
shown that they can solve certain problems much more 
efficiently than any classical device. Celebrated quantum 
algorithms have been built for the factorization of large 
numbers with exponential efficiency and for search 
in an unstructured database with a quadratic speed-up 
0. As first envisioned by Feynman in the 1980's, the 
simulation of complex quantum systems has also been 
shown to be more emcient on a quantum computer [jj. 

Here we study a quantum small-world network with 
disorder. We demonstrate the existence of a delocaliza- 
tion transition and investigate its dependence on disor- 
der strength, number of links and system size. We then 
build a quantum algorithm to simulate such a network on 
a quantum computer, and show that its efficiency signif- 
icantly overcomes classical computations. The algorithm 
is robust with respect to errors. 

We consider a circular graph with N = 2" r vertices. 
Each vèrtex is linked with its two nearest-neighbors. To 
this graph, pN shortcut links (connecting 2pN vertices) 
are added between random pairs of vertices (see an ex- 
ample in the inset of Fig.l) ||. A quantized version of 
this system with on-site disorder can be described by the 
N x N Hamiltonian matrix H = H + H x + H 2 . The 
first two terms give a one-dimensional tight-binding An- 
derson model well-known in sòlid state physics ||. The 



on-site disorder; Síj denote Kronecker symbols, and 
are independent random numbers whosc distribution is 
a Gaussian with zero mean and width W (the Gaus- 
sian is truncated at large vàlues). The matrix (-fíi)r, — 
V(Óíj+i + describes the links between nearest- 

neighbors, and (H 2 )ij = Efcíi Aj fc + s i,jk S 3,ik) 
the shortcuts which make the graph of small-world type, 
where {ik,jk} are the pairs of vertices connectcd by ran- 
dom links, and V = 1 is the hopping matrix element. 




diagonal matrix with entries (i?o)i 



describes 



FIG. 1: Level spacing statistics for H at n r — 14, p = 1/32, 
for three vàlues of the disorder: W = 0.5 (triangles), 1.3 
(empty squares) and 3 (full squares). The sòlid curves cor- 
respond to the Poisson distribution P(s) — e~ s and to the 
Wigner-Dyson distribution P(s) = (7rs/2)e _7rs ^ 4 . Number 
of disorder realizations (position of shortcut links and on-site 
disorder) is Nd = 10. Only the central half of the eigenval- 
ues is taken into account. Inset: a realization of small-world 
network with iV = 32 and p = 1/8. 

When p = 0, the system reduces to the one- 
dimensional Anderson model, for which all states are 
known to be localized. For small disorder, the local- 
ization length l varies as l oc 1/W 2 0- The addi- 
tional presence of shortcut links may induce delocaliza- 
tion. This can be checked through spectral statistics. 
Indeed, for localized systems, the eigenvalues are dis- 
tributed according to the Poisson distribution, provided 
the localization length is smaller than the system size. 
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FIG. 2: (color oniine) Evolution of the IPR £ with time t, 
for W = 0.5 and p = 1/32. Initial state is localized on one 
vèrtex. Curves correspond (from bottom to top) to n r = 
8, 10, 12. Each curve is shown as obtained by exact evolution 
(black lines), with Nd = 80, and by simulation by quantum 
gates (see text)(green/gray lines) with Nd = 100 and At = 
0.03. The arrows indicate the IPR at t — 2000 in the absence 
of shortcut links p — (exact evolution and simulation by 
quantum gates yield the same result, data not shown). 



for a wavefunction "F) = J2i It gives the number 

of vertices supporting the wavefunction (£ = 1 for a state 
localized on a single vèrtex, and £ = N for a state uni- 
formly spread over N vertices) . In Fig|2 and FigOU we 
display the time evolution of the IPR for a wave packet 
initially localized on one vèrtex. For W — 0.5, the sat- 
uration value grows with N in the presence of shortcut 
links, indicating that the wavefunction is no longer local- 
ized. On the contrary, for W — 3, the saturation value 
remains close to its value in the absence of links and does 
not change significantly with N, implying that the sys- 
tem is still localized. In a more quantitative way, Fig0] 
presents the saturation value of the IPR as a function 
of n r for different vàlues of p. The data confirm that at 
W = 3 the system remains localized. On the contrary, 
a clear delocalization is visible in the presence of short- 
cut links for W = 0.5. The data are in good agreement 
with the law £ cx N a , with a « 0.58 for p = 1/32 and 
q « 0.84 for p = 1/16 (the maximal value a = 1 is ob- 
tained at p — 1/2, data not shown). This shows that the 
delocalization transition for p = 1/16 and p — 1/32 takes 
place approximately at W ~ 1. In the limit of weak dis- 
order Wcl, the transition is expected to take place at 
smaller vàlues of p oc W 2 [3j . 
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FIG. 3: (color oniine) Same as FigUJwith W = 3. 



On the contrary, in the delocalized phase the eigenval- 
ues follow the Wigner-Dyson distribution corresponding 
to Random Matrix Theory, which generally characterizes 
quantum chaotic systems and ergodic wavefunctions @- 
Our numerical diagonalization of H at fixed p shows a 
transition from Poisson to Wigner distribution as W de- 
creases. A typical example is shown in Fig.l at p = 1/32 
and W = 3 (localized phase), W = 1.3 (intermediate 
statistics), W = 0.5 (delocalized phase). This indicates 
that a delocalization transition takes place in this system. 

The localization properties of this quantum system can 
be analyzed more precisely through the Inverse Partici- 
pation Ratio (IPR), defined by £ = |* i | 2 /E l l**| 4 




FIG. 4: IPR at time í = 2000 as a function of n r for W = 0.5 
(full symbols) and W = 3 (empty symbols), for p — 1/16 
(triangles), p — 1/32 (squares) and p — (circles), with 20 < 
Nd < 160. Initial state is localized on one vèrtex. Straight 



lines correspond from top to bottom to £ = 0.78ÍV 



£ = 



3.13iV u 



£ = 70.15, £ = 13.55, £ = 7.91, £ = 5.69, with 



N = 2™ r . Logarithm is decimal. 

This system can be simulated on a quantum com- 
puter, using 0((log N) 2 ) quantum gates for a network of 
N = 2" r vertices, and n q = 3rt r + 3 qubits. We start from 
an initial wave packet encoded on the quantum registers. 
For example, the initial one-vertex states used in FigsJJJ- 
Ulcan be constructed efficiently from a state localized in 
the ground state of the quantum computer by at most 
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n r single-qubit flips. Our quantum algorithm performs 
the evolution of the wave packet by slicing the prop- 
agator exp(ïifí), using the relation e l ( H a+Hi+H 2 )&t _ 
e iH ^ e iH^ e iH % ò.t e iH^ e iH ^ + (At 3 ) for a short 
period of time Aí (see e.g. 0, E3). Each unitary op- 
erator is then simulated by quantum gates. We use 
in particular rotations on the j-th qubit by an angle 
(f>/2: Rj((f>) = exp(i(/>cr|/2) (a z being a Pauli matrix); 
controlled-not operations CNOTjj, that is bit-flip on the 
j-th qubit conditioned by the i-th qubit; multi-controlled 
rotations Cii,...,i tl ,rn,...,ri li ,j(0)i that is rotations by an an- 
gle 8 on the j-th qubit if and only if the qubits ik takes 
the value r\k <E {0, 1} for 1 < k < \i. 

• The transformation \ï) — > e lH ° At \i) consists in mul- 
tiplying each basis state \i) by a Gaussian random phase 

exp(ieiAí). For some integer n s , and a = WAt^J~- ^ , 

let us choose randomly n r + n s angles 4>k, 1 < k < n r , 
and <f>'u, 1 < k < n s , independent and uniformly 
distributed in [—a/2, a/2]. Each e^Aí is replaced by a 
random variable ±0i ± (f>% ± • •• ± (/)'s-i ^ 0si which for 
large n s tends to a Gaussian random variable of width 
WAt. This can be simulated by applying the operator 

for some value of n e . The ik and jk are chosen randomly 
between and n r — 1. This step requires (3n s + n r ) 
gates. 

• To perform the transformation \i) — > e líílAt |z), 
we first apply a quantum Fourier transform (QFT) 
to turn it into the diagonal transformation |fc) — » 
exp (2iAí cos 2^) |/c). Following ^(Ji we introduce the 
operator iï 7 ((9) = HS l He~ l i< HS^He^ï^ HS 1 with 
£" m = üj=2 ( 7rm /2"'~ 1 )- can be shown that 
e -; 7 cose = R i/2 ((j)R l/2 (-6) + 0( 7 3 ) for small 7, with 
8 = 8 — nai if #/27r = 0. 010,2 . . . a n . The diagonal oper- 
ator can thus be approximated by exp (2iAícos^£) = 

{R 1 / 2 {8)R 1 / 2 (~8)) L + 0(Lj 3 ), with L an integer and 7 a 
small parameter chosen such that = — 2iAt. We then 
perform an inverse QFT. The two QFT require n r (n r + l) 
gates, and the simulation of the diagonal term requires 
2L(5 + 3n r ) gates. 

• The transformation \i) — > e tH2At \i) acts on the sub- 
space spanned by \ik) and \ jk), where ik and jk are linkcd 
by a shortcut link, through the 2x2 submatrix e lAta . 
Let us first assume that p is of the form p = 1/2 P . For 
fx = p — 1, the operator C^^..^ Vlj ... tV j(8) acts on the 
2™ r_A1 basis vectors whose qubits ik, 1 < k < fi, are re- 
spectively equal to 77^-: it corresponds to the creation of 
2»i,.-p-i = 2" r ~ p = pN links. In order to have less regu- 
lar shortcut links, we first perform a permutation on the 
vertices. To do this, we randomly choose n p integers Ofc 
and bk, for some integer n p . It is better to take the ak 
in [0.2ÍV, 0.8-ZV] and odd. Then we define the operators 
Uk\ï) = |(afeï + bk) mod N) and the inverse operators 
Vk\i) = \o-k 1 (i — bk) mod N). A permutation can be sim- 
ulated by the sequence of gates P = Y\.kLi Uk CNOT ifc j h 



where the ik and jk are chosen randomly. Application of 
the permutation P, followed by a multi-controlled rota- 
tion C , ill ... )<MjlJll ... )W (At) and P~\ gives é H * At . The i k 
and rjk in the controlled rotation are also chosen ran- 
domly. In the general case, where p ^ 1/2 P , we ex- 
pand pN in base 2, such that pN = ^2 Pfc . Then we 
replace the multi-controlled gate in the above descrip- 
tion by a multi-controlled gate for each pk appearing in 
the decomposition of pN. This gives a sequence of gates 
C>) A") „(*) _(*) ,■(«,) (Aí), where fj, k = n r - p k - 1, 

and the éy 1 , jW and rjíy are chosen randomly. Each 
operator U k consists of a multiplication and an addition 
modulo N, which can be performed using (2n r + 3) an- 
cilla qubits and O(n^) quantum gates [llj. Each multi- 
controlled gate can be perf ormed by O(n^) Toffoli, Cnot 
and single qubit gates |L2|. 

In total, the simulation of a network of N = 2™ r ver- 
tices for one unit of time with fixed parameters Aí, n s , 
L and n p can be done by this method with 0(n^) quan- 
tum operations and 3n r + 3 qubits. Classically, a similar 
method can only be implcmcntcd in O(N) operations at 
best. The quantum simulation is therefore exponentially 
faster. This remains the case even if the parameters n s 
and n p arc allowcd to grow lincarly with n r to improve 
accuracy (the cost becomes 0{n 3 ) quantum gates). 

The algorithm simulates the small-world network efh- 
cicntly but at the cost of several approximations. In order 
to check its convergence and accuracy, we implemented 
it on a (classical) computer. In Figs l2l3l we display the 
result of this computation for the parameters Aí = 0.03, 
n s = 30n r , L = 10, and n p = 3n r alongside the exact 
evolution, showing that the algorithm is quite accurate 
for these vàlues, and enables to monitor precisely the 
delocalization transition with good accuracy. The com- 
putation accuracy is not very sensitive to fixed vàlues of 
L and Aí: the total size N can be changed by orders of 
magnitude (factor of 64 in our case) without modification 
of thcsc parameters. 

To estimate the total complexity of the algorithm, we 
should take into account the number of quantum mea- 
surements and the number of iterations of the map. In 
order to sec the delocalization transition, it is sufncicnt 
to estimate the spreading of the wavefunction, which can 
be done by a constant number of quantum measurements 
[ï^ . Still, the initial wave packet should have enough 
time to spread in order for the localization length to be 
estimated. For the parameters of Fig[21 we determined 
the time r needed for the IPR to reach half of its max- 
imal value. In the delocalized phase for p < 1/2, our 
numerical results give the sealing t oc A' 3 with j3 « 0.83 
(p = 1/16) and « 0.69 (p = 1/32) (data not shown). 
This mcans that the total cost of the quantum algorithm 
will scale as 0(N 13 ), compared to 0(N@ +1 ) for the clas- 
sical one (dropping logarithmic factors). This implies a 
better than quadratic gain for the quantum computation, 
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but no exponential gain. In contrast, for 1/2 < p < 2 we 
find that r w logAT (data not shown) 0|. In this case, 
the algorithm may reach exponential efíiciency and en- 
able to perform precise studies of this percolation-likc 
transition for very large vàlues of N. The exact algo- 
rithm complexity depends on the properties of the phase 
transition near critical W value. 
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we took the part of the algorithm which generates the 
random shortcut links as exact, all other parts being per- 
formed with errors. The results displayed in FigElshow 
that with moderate levels of imperfections (e rí 10~ 7 ) 
the simulation of the small-world network is very close 
to the exact computation, in absence of any quantum 
error correction. 

In conclusion, we have shown that quantum disor- 
dered small-world networks, which display a delocal- 
ization transition, can be simulated more efficiently on 
quantum computers than on classical ones. The al- 
gorithm can be performed accurately on realistic few- 
qubit quantum computers in presence of moderate error 
strcngth. 

We thank A. Pomeransky and O. Zhirov for helpful 
discussions. We thank the IDRIS in Orsay and CalMiP in 
Toulouse for access to their supercomputers. This work 
was supported in part by the project EDIQIP of the IST- 
FET program of the EC. 
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FIG. 5: Evolution of the IPR with time in the presence of 
static errors, for n r = 10, W — 3 and p = 1/32, with Nd 
I Üí ! (static errors are the same for all realizations) . Initial 
state is localized on one vèrtex, Aí = 0.03. From top to 
bottom: e = 10" 4 , 10" 5 , 10" 6 , 10" 7 and 0. The arrows mark 
Ç at í = 500 for W — 0.5 and same other parameters, with 
from top to bottom e = 10" 5 , 10" 4 , 10" 6 , 10" 7 and 0. Inset 
shows the three lowest curves on a different scale. Data from 
e = 10 -7 are indistinguishable from e = 0. 

These results show that a perfect quantum computer 
gives a significant gain in the simulation of quantum 
small-world networks. However, realistic quantum com- 
puters are prone to errors and imperfections. It is there- 
fore important to test the resilience of the algorithm to 
such effects. In FigGD we show the result of numerical 
simulations of the algorithm in presence of errors. The 
error model chosen corresponds to static imperfections. 
These errors can exist independently of the coupling with 
the external world, and have parametrically larger effects 
than random noise in the gates ^J. Between each gate 
the system evolves through the additional Hamiltonian 
He = J2i ^i°í +2j Ji a ï a i+i > where the second sum runs 
over nearest-neighbor qubit pairs on a circular chain. The 
Si are randomly and uniformly distributed in the interval 
[-S/2, S/2]. The couplings Ji represent the residual static 
interaction between qubits and are chosen randomly and 
uniformly in the interval [— J, J]. We suppose that each 
gate in the quantum algorithm is instantaneous and sep- 
arated by a time T g during which He acts. We take 
one single rescaled parameter e which describes the am- 
plitude of these static errors, with e = 5r g — Jr g . In 
the numerical simulations, to save computational time 
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